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SUMMARY 

When adapting a vehicle for use by a handicapped driver, it is often 
necessary to employ a high gain steering controller to compensate for limit- 
ations in the driver's range of motion. Because such a driver /vehicle system 
can become unstable as vehicle speed is increased, it is desirable to use a 
computer simulation of the driver /vehicle combination as a design tool to 
investigate the system response prior to construction of a controller and 
road testing. While there are a number of different driver models in exis- 
tence, they all contain some unknown driver parameters which must be iden- 
tified prior to use of the model for system analysis. This work addresses 
a means to collect the data necessary for identification of these driver 
model parameters without extensive instrumentation of a vehicle to measure 
and record vehicle states. 

The procedure consists of three steps. First, a road test is conducted 
with the driver in a normal vehicle, during which only the steering wheel 
angle and the vehicle speed is recorded. Next, the data from the road test 
is input into a computer n. ^del of the vehicle which - ’crates the vehicle 
equations of motion with the given speed and steering 4.. its to yield the 
vehicle states, some of which the driver senses. Finally, with the sensed 
vehicle states as inputs and the recorded steering response as output, a 
least squares parameter identification procedure is used to coi “'ute the 
parameters in the proposed driver model. 

Initial tests of the procedure identified all of the driver parameters 
with errors of 6% or less. 
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IMTRODUCTION 


The Rehabilitation Engineering Program at Texas A&M University is cur- 
rently involved in the evaluation and design of low effort - high gain auto- 
motive control devices for handicapped drivers. Experience has shown that 
some vehicles with high gain steering are relatively easy to control at mod- 
erate speeds and others require maximum driver effort to maintain control at 
very low speeds. In order to better understand the causes of this wide vari- 
ation in handling proterties and to quantify the effects of changes in various 
steering system parameters on vehicle response, development of a computer sim- 
ulation of the driver-vehicle combination was begun. 

At present, there are a number of vehicle models [3] and driver models 
[1,5,6] that are available. Use of a typical vehicle simulation requires know- 
ledge of the geometry, inertial properties, and tire characteristics of the 
vehicle to be studied. Most of these can be obtained by direct measurement 
or are easily estimated. 

As with the vehicle models, most of the driver models were found to have 
several coefficients whose magnitude is dependent on the characteristics of 
the driver or his environment. In order to use a driver model, a means to 
quickly and inexpensively identify these unknown parameters in a driver model 
is required. Classically, this identification is accomplished by running a 
road test and recording the driver inputs and all of the vehicle motion vari- 
ables that the driver may respond to. As illustrated in Figure la, the 
vehicle motion data is used as an input for the driver model and the driver 
response data is compared to the results of the driver model to generate an 
error. There are a number of parameter identification techniques available 
[4] to determine the coefficients in the driver model that will minimize the 
error in predicted and measured driver response. The difficulty in employing 
this procedure lies in the instrumentation required to record the vehicle 
motion variables. Variables such as heading angle or yaw rate require a 
gyroscopic device which is expensive and bulky. The position of the vehicle 
on the road may be obtained with an optical tracking device or by integrating 
the output of accelerometers on an inertial platform. Either of these methods 
is expensive and the equipment is not easily moved from one vehicle to another. 
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Figure la. Typical Parameter Identification 
Procedure 
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Figure lb. Proposed Procedure 
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In order to minimize the instrumentation necessary to obtain the required 
data, an alternate procedure was examined. As illustrated in Figure lb, the 
only measured variables are driver inputs (steering wheel angle and vehicle 
forward velocity). These inputs are then used with a vehicle simulation to 
determine what the motion must have been during the test. In this manner the 
required data for driver model identification is obtained with a minimum of 
instrumentation. 

VEHICLE MODEL 


The vehicle model used in this procedure was chosen because it was con- 
sidered to be the simplest model available that exhibited the handling pro- 
perties under study. Figure 2a is an illustration of the vehicle model. This 
vehicle has three degrees of freedom including the forward and lateral posi- 
tion of the mass center P, (r and r ) and the heading angle (i()) . Two 

r A I I 

coordinate systems are employed. The X, Y, Z system is fixed to the roadway 
and has the associated unit vectors I, J, and K. The x, y, 2 system is fixed 
to the car and has unit vectors i, j, and k. 

Using these definitions, the acceleration of the mass center, P, may be 
shown to be 

COS + tpY Sin i + [rp^ cos - ilp^^ sin 4>] j (1) 

The free body diagram for the vehicle is illustrated in Figure 2b. In 
order to simplify the analysis, secondary forces such as tire rolling resis- 
tance, self aligning torque, aerodynamic drag, and gyroscopic moment are con- 
sidered negligible. Applying Newton's Laws to the free body diagram results 
in the following equations of motion. 

*L * ’’fR ‘r * ''rl * ’’rr ■ ”'py ♦ - "Vx 

\ ^fR ^R^ ~ 

+ (a^ - a^) sin 6^ - (3) 
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A third equation may be obtained by recalling that the vehicle forward 
velocity is known. The x component of a^, therefore, may be expressed as 

ifpjj cos + i-pY sin ij) « V (A) 

Equations (2), (3), and (4) represent the three equations of motion for 
the vehicle. Unfortunately, the tire forces are unknovm 

and must be related to the vehicle motion parameters. It is well documented 
[2] that a rolling pneumatic tire under the Influence of a lateral load 
exhibits a viscoelastic deformation of the tread surface. This deformation 
results in a deviation of the wheel center velocity from its expected direc- 
tion. This deviation is termed the slip angle and may be expressed in terms 
of the vehicle motion variables. Using the defined coordinates, Tpy, 

and i)») the slip angle at each wheel may be shown to be: 


fL 


6 - tan 

Lt 


*fL 


RL 


= 6 . 


-r- 

-f 


cos - r„v sin ’I' + Sj’l' 


PX 


V + (aj + a^)i|> 

• • 

COS 4> - sin + aj\(/ 


PX 


= - tan 




V - (a, - ajil. 


PY cos ij; - tpj^ sin - a^ilil 


RR 


tan 


V + (aj + a^)4» J 

cos - Tpjj sin 4» - a^4>1 




V - (a^ - a^)^» 


(5) 


( 6 ) 


(7) 


( 8 ) 


where 6 and 6 are left and right steering angles as defined in Figure 2a. 
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Once the slip angle at a given wheel is specified, the tire lateral force 
may be determined. For this work, the tire used is a Goodyear FR70-14. 

Figure 3 illustrates the relationship between slip angle and lateral force for 
various normal loads on the tread surface. 

Given the equations of motion (2-4), the slip angle e’^^pressions (5-8) 
and the tire force function illustrated in Figure 3, the motion of the vehicle 
may be determined if the forward velocity V, and the front wheel angles, 6^ 

Assuming Ackermann Steering, 6^ and 6^ may be related to 


and 6 are known 

K 

the steering wheel angle by the following relationships 

L 

X 


tan (Kg^6) 


(9) 
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6 = tan ^ 7 - ^ 

R (x - a^) 


6, = tan“^ 7 — V- 

L (x - aj) 


( 11 ) 


where represents the steering gain, 6 is steering wheel angle, and L is 
wheelbase. 

Now, given the driver inputs of steering wheel angle, 6 , and forward 
velocity, V, the motion of the vehicle may be determined. Because of the non- 
linearity of th .otion and constraint equations, the most expedient means of 
obtaining a solution is by utilizing a numerical integration procedure on a 
digital computer. The computer that was used in the following example was a 
PRIME 730 minicomputer. The program was written in FORTRAN and utilized a 
predictor>-corrector numerical Integration scheme. 


VALIDATION OF VEHICLE MODEL 


In order to verify the accuracy of the vehicle simulation, a series of 
test runs was made with a vehicle that was instrumented to record several of 
the vehicle motion parameters as well as the driver responses. The driver 
responses were used as input data for the vehicle simulation The results of 
the vehicle simulation were then compared to the vehicle motions recorded 
during the test. Several test runs were made at different vehicle speeds and 
over different courses. Figure 4 illustrates the comparison between predicted 
and measured heading angle, yaw rate, and lateral acceleration as the vehicle 
passed through an offset alley. As shown in Figure 4, the vehicle simulation 
provides a good estimation of the vehicle motion. The large difference in 
measured and predicted heading angle is due to a difference in reference 
position. The simulation automatically sets the vehicle’s initial heading to 
zero degrees while the measurt^d heading depends on the vehicle orientation 
when the gyroscope unit is switched on. When this bias is removed, the 
results compare quite well. 
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Figure 4. Comparison of Measured and 
Simulated Vehicle Motion 



VALIDATION OF IDENTIFICATION PROCEDURE 


The driver model used to validate the parameter identification procedure 
assumes that the driver can be represented as a two part cascade of a brain 
sponse and 4 . neuromuscular lag. Reaction time delay and precognition or 
preview are not Included since they tend to cancel and, at any rate, the par- 
ticular form of the driver model is not critical to the identification pro- 
('(•d’lre at this stage. The inputs to the driver model are the lane position 
irror, and the handing angle error, . For a straight roadway = r^^ 

and e = The output of the brain is the commanded steering wheel angle 6^ 


and is modeled as 


K (Tpy + nO 


( 12 ) 


where K is the brain gain and n is a weighting factor. The brain output is 
the input to the neuromuscular system which is modeled as 


6 = ^ (6' - 6) (13) 

where t is the neuromuscular time constant. Thus, the parameter Identifi- 
cation procedure must compute values of n, and t for which the computer 
model output best fits the measured data in a least sijuares sense. 

Since the driver parameters K, and r are not available from an actual 
dr iver/vGhlcle test a computer simulation was used to generate test data that 
allows direct comparison of computed parameter values to true parameter 
values. The test data was generated using the same driver model as the 
identification code, but a different vehicle model. The vehicle model used 
was more complete a.id had been thoroughly checked out. 

The test maneuver was lane tracking on a straight road at a nominal speed 
or 55 mph with an initial lane error of five feet. As shown in Tabl^ 1 the 
largest error between the true and computed parameter values is 6%. 


Parameter 

True Value 

identified Value 

Percent Error 

X 

0.1 

0.100 

02 

K 

- 0.1 

-0.094 

-62 

n 

100 

94.40 

-5.62 


Table 1 

Errors in Identified 
of the Driver Model 

Values 

Parameters 
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CONCL\ SIGNS 


As sbovm in Figure 4, the computer model provides a very accurate simu- 
lation of the actual vehicle states. This is a critical result, since the 
data collection scheme proposed in this paper is based on the premise that 
vehicle states, which are inputs to the driver, can be accurately constructed 
from simulation, thuc eliminating the need to instrument each vehicle to oe 
tested with an inertial platform. Also, assuming that an appropriate driver 
model is used, the identification procedure accurately computes the driver 
model parameters as indicated by the data in Table 1. Thus, each part of 
whe overall procedure has been independently validated. 

The final test of th . procedure will be identification of driver model 
parameter's from an actual road test and comparison of computed vehicle states 
and driver outputs to actual measured v lues. This testing is currently in 
progress. 
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